INF-SUP ESTIMATES FOR THE STOKES PROBLEM IN A 
PERIODIC CHANNEL 

JON WILKENING * 

Abstract. We derive estimates of the Babuska-Brezzi inf-sup constant j3 for two-dimensional 
incompressible flow in a periodic channel with one flat boundary and the other given by a periodic, 
Lipschitz continuous function h. If h is a constant function (so the domain is rectangular), we show 
that periodicity in one direction but not the other leads to an interesting connection between f3 
and the unitary operator mapping the Fourier sine coefficients of a function to its Fourier cosine 
coefficients. We exploit this connection to determine the dependence of /3 on the aspect ratio of the 
rectangle. We then show how to transfer this result to the case that h is C 1,1 or even C ' 1 by a 
change of variables. We avoid non-constructive theorems of functional analysis in order to explicitly 
exhibit the dependence of /3 on features of the geometry such as the aspect ratio, the maximum slope, 
and the minimum gap thickness (if h passes near the substrate). We give an example to show that 
our estimates are optimal in their dependence on the minimum gap thickness in the C case, and 
nearly optimal in the Lipschitz case. 
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1. Introduction. Many problems of industrial and biological importance in- 
volve fluid flow in narrow channels with moving boundaries [7J 111) . Examples include 
the flow of oil in journal bearings or between moving machine parts, the flow of air be- 
tween disk drive platters and read- write heads, or the flow of mucus under a crawling 
gastropod [T3J . A primary objective in all these problems is to solve for the pressure 
required to maintain incompressibility. Indeed, it is the pressure that determines the 
load sustainable by a journal bearing, and that provides propulsion against viscous 
drag forces in peristaltic locomotion. However, only the gradient of pressure enters 
directly into the Stokes or Navier-Stokes equations; thus, regardless of the method 
used to solve the equations, the pressure must be determined via its gradient. 

The fundamental fact that makes it possible to extract p from Vp is that the 
gradient is an isomorphism from L^(fl), the space of mean-zero square integrablc 
functions, onto the subspace of linear functionals in i/ _1 (J7) 2 that annihilate the di- 
vergence free vector fields u e H^tt) 2 ; see Section [5] below. The inf-sup constant 
/3 (or rather, its inverse) gives a bound on the norm of the inverse of this operator. 
Thus the magnitude of p (and our ability to estimate errors in p) depends to a large 
extent on the size of However, to the author's knowledge, every existing proof 
(e.g. [HIE]) that is finite relies on Rellich's compactness theorem to extract a sub- 
sequence whose lower order derivatives converge, making it impossible to determine 
how large /3 _1 might be or how it depends on fi. The proof in [4] also uses the closed 
graph theorem, which, like Rellich's theorem, leads to constants that depend on SI in 
an uncontrollable way. These proofs are appropriate for pathological domains with 
bulbous regions connected by thin, circuitous pathways; however, for "nice domains" , 
it should be possible to obtain better estimates of the constants — existing theorems 
are of limited practical use. 
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Au + Vp = f 
V • u = 



h{x) 



max /i(a; 

0<x<L 



h 



min /i(x) > 

0<x<L 



u = 



x = 



x = L 



Fig. 2.1: Two dimensional Stokes flow in a periodic channel. The left and right 
boundaries have been identified and are considered to be part of the interior of the 
domain. 



In this paper, we derive explicit estimates of the inf-sup constant /3 for two- 
dimensional incompressible flow in a periodic channel with one flat boundary and the 
other given by a periodic, Lipschitz continuous function h(x). Our goal is to determine 
how /3 _1 depends on features of the geometry such as the aspect ratio, the maximum 
slope, and the minimum gap thickness (if h passes near the substrate). Although 
these requirements on fi are fairly restrictive, such geometries do cover a wide range 
of interesting applications. 

Our interest in this problem arose in the course of deriving a-priori error estimates 
for Reynolds' lubrication approximation (and its higher order corrections) with con- 
stants that depend on fl in an explicit, intuitive way; see [H] and also J7HTU1III] for 
background on lubrication theory. These a-priori estimates were used by the author 
and A. E. Hosoi to monitor errors in the lubrication approximation while studying 
shape optimization of swimming sheets over thin liquid films; see |13) . 

2. Preliminaries. In this section we briefly review the weak formulation of the 
Stokes equations, emphasizing the role played by the Babuska-Brezzi inf-sup condi- 
tion; see e.g. [21 El [6] for a more detailed account. 

Consider the two-dimensional, cc-periodic Lipschitz domain Q shown in FigureJ^HJ 



il = {{x,y) : x£T, 0<y<h(x)} 7 heC a ^(T), T=[0,L] p . (2.1) 



The case of non-zero Dirichlet boundary conditions may be reduced to the homoge- 
neous case by subtracting off an appropriate function to transfer the inhomogeneity 
from the boundary conditions to the body force f; see e.g. [3]. We treat f2 and T as 
C°° manifolds by identifying the points 



and adding a coordinate chart to each that "wraps around" . In particular: a function 
in C k (Q) or C k (T) is understood to have k continuous periodic derivatives; 9f2 = 
To U Ti; dT = 0; the support of a function (j> E C^(S1) vanishes near To and Ti but 
not necessarily at x = and x = L; and the Sobolev spaces H k (n) and Hq(£1) are the 
completions of C k (Q) and (f2) in the || • \\k norm, and thus contain only ^-periodic 
functions with appropriate smoothness at x = 0, L. 

In the weak formulation of the Stokes equations, we seek the velocity u and 



n : (0, i/) ~ 
T : - L 



0<y<h(0), 



(2.2) 
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pressure p in the spaces 

X = Hl{tt) 2 1 M = Lg(£l) = jp e L 2 (Q) : J pdA = oj, (2.3) 



respectively, such that 



a(u, v) + 6(v, p) = (f , v) (2.4a) 
b(u,q) =0 (2.4b) 



for all v € X and q e M, where the body force f may be any linear functional in the 
dual space X' — 1 (SI) 2 and 

a(u,v)= / Vu:VveL4, b(u,p) = - pV-udA. (2.5) 
,/n is! 

We endow M with the L 2 norm || • ||q and X with the energy norm (i.e. the H 1 semi- 
norm) ||u|| a = y/ a(u, u), which is equivalent to the H 1 norm ||u||i = vIMIo + ll u llo 
due to the Poincare-Friedrichs inequality (see Lemma IA.2|) : 

||u|| < -i||u|| , (ueX), hi=maxh(x). (2.6) 

V8 0<a;<L 

Next we define the operators B : X — » M' and B' : M —> X' via 

(Bu,p) = 6(u,p) = (£'p,u), (£ = div, B'^grad). (2.7) 
B and B' are clearly bounded and satisfy 

||B|| = = sup sup |b [ U '^| < s/2, (2.8) 

where M = M \ {0} and X = X \ {0}. We note that ifuel then V • u € M, i.e. the 
divergence of u has zero mean; hence, 

V :=\ictB = {ue X : V ■ u = 0}. (2.9) 

The Babuska-Brczzi inf-sup condition 

3/3>0 such that inf sup J^"',^ > (2.10) 
P eM nex IbllolMla 

is precisely the condition required for B' to be an isomorphism onto its range with 
inverse bounded by || (_B') _1 || < /3 _1 . Once we know the range of B' is closed, we may 
take the polar of the equation ran(_B')° = ker(B) = V to conclude 

ran(B') = V° = {t £ X' : (f , u) = whenever u G V}. (2.11) 

As ^° is naturally isomorphic to (X/V)', we see that B : X/V -> M' :(u + V)^Bu 

is the adjoint of the composite map M -—> V° (X/V)', and is therefore itself an 
isomorphism with the same bound on the inverse. Identifying X/V with 

^^{ueX : a(u, v) = whenever v e 7}, (2.12) 
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we learn that the restriction of B to V- 1 is an isomorphism onto M' , which would 
be essential to the analysis of the Stokes equations if the right hand side of ()2.4b|) 
were inhomogeneous. Other interesting solutions of Bu = if with (p G M' (requiring 
e.g. u e L°°(S1) 2 n X or V x u = rather than u G V^-) are studied in [T|. Finally, 
we define A : X -> X' and I : V -> V via 

(Au, v) = o(u, v), (u, v G X), (Au, v) = a(u, v), (u, v G V). (2.13) 

Both are isometric isomorphisms in the |j • ||„ norm. 

The weak solution (u,p) of (|2.4[) must satisfy u e 7 so that Bu = 0. But then 
Au + B'p = f requires 

(*) Au = f, (f) B'p = f-Au, (2.14) 

where f = f|y £ V and we note that (f — Au) G range(B') = V° iff u satisfies (*). 
Since A and 2?' are isomorphisms onto their ranges, a unique solution of (|2.4p exists 
and we have the estimates 

||u||„ = \\f\\ v > < \\f\\x> = sup Ij^l, ||p|| <2/3- 1 ||f||^. (2.15) 

In summary, the inf-sup condition (|2.10p is the key to analyzing the weak formulation 
of the Stokes equations — it is equivalent to the assertion that the gradient B' is an 
isomorphism from M = L 2 ^ (f2) onto the polar set V° of linear functionals in X' that 
annihilate the divergence free vector fields u E V. 

It is instructive to compare the inf-sup condition written in the form 

/?IW|o < \\B' P \\ X ' = ||Vj»||_i < V2\\p\\ (pe£|(n)), (2.16) 
to the Poincare-Friedrichs inequality for mean-zero functions: 

Ibllo < C||Vp||o => (l + ^r^llplli^llVpllo^lNIi (pe4(fi)). (2.17) 

Whereas (|2.17[) is easy to prove for p G Hq(SI) (with C = in our case), it is 

more challenging to prove for mean zero functions p G H^(£l). The usual proof [2j[5] 
relies on Rellich's theorem that i? 1 (S7) is compactly embedded in L 2 (Q). As a result, 
the proof does not tell us how large the constant C might be or how it depends on fi. 
Similarly, the usual proof [4] of (|2.16j) makes use of Rellich's theorem that L 2 (fl) is 
compactly embedded in iJ _1 (f2) = Hq(CI)'; however, there is an added complication 
not present in proving (|2.17[) : it must first be established that 

Mo < CdMU + UVplU), (peI 2 (S])). (2.18) 

This can be done in our case (if h G C 1:1 (T)) by flattening out the boundary and 
constructing appropriate extension operators from £f -1 (f2) to if _1 (T x M) to reduce 
the problem to a case that can be solved using the Fourier transform; see Duvaut and 
Lions [J and also Nitsche [3], who used a similar technique to prove Korn's inequality. 
In this paper, we show how to bypass (|2.18[) and prove (|2.16[) directly without invoking 
Rellich's theorem, which allows us to determine how the constant (3 depends on f2. We 
present two versions of the proof: one assuming h G C ' (T), and the other assuming 
only that h G C 0,1 (T), i.e. that h is a periodic, Lipschitz continuous function. Our 
proof does rely on the boundary of being the graph of a function h(x); however, we 
feel this is a sufficiently important case to warrant a separate analysis. We sketch a 
proof of (|2.17p that avoids Rellich's theorem in Appendix [Bl for comparison. 
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3. A rectangular channel. In the following theorem, we prove that B' in (|2.7p 
is an isomorphism onto its range (with (3 = |min(l,4-^)) when Q is the x-periodic 
rectangle R = T x (0, H) of height /f. In Sections 0] and 03 we will transfer this 
result to a general x-periodic domain £1 by a change of variables. It is useful in this 
change of variables to know that the constant Ci in Theorem 13. II (and especially in 
Corollary OQJ does not diverge as H approaches zero. 

The periodicity of the domain in one direction but not the other leads to an 
interesting relationship between the inf-sup condition and the unitary operator map- 
ping the Fourier sine coefficients of a function of one variable to its Fourier cosine 
coefficients. By studying this operator, we can obtain explicit estimates of (3 and its 
dependence on L/H. 

Recall that every u E H^(R) must be zero (in the trace sense) on the top and 
bottom walls but not necessarily on the side walls, where it is only required to be 
periodic. Such a function can be expanded in a sine or cosine series in the y-direction 
and differentiated term by term. (If u € H 1 (R) is not zero on the top and bottom 
walls, only the cosine series can be differentiated term by term). 



Theorem 3.1. For all q e Lg(R), 



kilo ^CillcUH 2 .! +C7 2 R<?H 



(3.1) 



where C\ = max f 9, ^jp), C 2 = 9, and ||/||_i =sup„ 6H i (K) 



Proof. We may expand any q 6 Lq(R) and u <E Hq(R) in a Fourier series 
( a n j V2 



q(x,y) 


-£ 

n£Z 




u(x,y) 


-£( 

ngZ 


C„0 " 


so that 






nil 


= E 

ZxN 


LH\ 


\\< 


= E 


LH 



00 



3=1 



7r jy\ 

{0S ir)' 



2-ninx 



E 

n£Z 



3=1 



2-ninx 



3=1 



'raj I 



3 nj V2cos- l 



: < 



E 

raSZ 



\/2 sin 

3=1 



IT/ 



/2nn\ 2 /7T7 

(— ) + (i 



-"3 I 



= E LH \ b 

ZxN 

= E^ 



(3.2) 



27m \ 2 



Here No = {0}UN and the sums are over ordered pairs (n, j). 
Z x N \ {(0, 0)}. We claim that 



E 

(ZxNo 



E 



i 

LH{2im/L) 2 \a nj \ 2 
i (2Trn/L) 2 + (TTj/H) 2 

VI 

LH(2irn/L) 2 \b nj \ 2 
(2nn/L) 2 + (itj/H) 2 



\d x q\\ 



(3.3) 

Let us denote (ZxN )' = 



E 

ZxN 



E 

ZxN 



LH(jrj / H) 



{2-kti/L) 2 + (nj/H) 2 
Al 



(3.4) 



LH(iTj/H) 2 \h 



By 



(27TO/L) 2 - 

T 

B 2 
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Here A\, A 2 , Bi and B 2 are labels to represent the indicated sums. The horizontal 
assertions clearly hold (since q G L\ (R) =>■ aoo = 0) while the vertical assertions follow 
from the Cauchy-Schwarz inequality and a particular choice of u to show that two of 
the upper bounds are least upper bounds: 



(d x q,u) = ( (q)(-d x ii) dA = ^ LHb- nj — —j—d nj ^j < B\ /2 \\u\\ a 

•' R ZxN 

(d x q,u) = ^ £# a -«j ( - — ^ c "j) < ^i /2 ||"l|a 

(ZxN )' 

(d y q,u) = / (q)(-d y u)dA = ^ LHa_„j(- ~^dnj) < A^WuWa 

J R ZxN 

ZxN 

{2mn/L)b- n j 1/2 ,„ , 

d n j = /r ^ 2 : 73; 7%^ IKIU = - B i ? {d x q,u)=Bi, (3.5) 



(27TO/L) 2 + (ttj/H) 2 

-(ttj/H) a_ niJ - ^ |Uf|| _ 4 i /2 
{2-Kn/Lf + (ttj/H) 2 



^= r,Jn^' t - n ^ => HI- = 4", (^,«>=^. (3-6) 



The choices of c ra j analogous to (|3 . 5[) and (|3.6[) do not generally lead to functions u 
that satisfy the boundary conditions on the top and bottom walls; hence, we cannot 
replace the inequalities in ()3.4j) by equalities. 

The theorem will be proved if we can show that 

6(Ai + A 2 ) + (1 - 6){Bi + B 2 ) < C X B X + C 2 A 2 (3.7) 

for some 9 £ [0, 1]. The result (|3 . 15[) below turns out to be independent of 8, so we 
set 9 = 1 here for simplicity. We will prove (|3. T[) by slicing the lattices (Z x No)' and 
ZxN into vertical strips and showing that 

A lin < dB lin + (C 2 - l)A 2 , n , (n £ Z), (3.8) 

where the subscript n indicates that only the terms in strip n should be included in 
the sum, e.g. A 1>3 = £~ LH(6ir/L) 2 \a 3j \ 2 /[(6ir/L) 2 + (ttj/H) 2 }. Since A 1>0 = 0, 
the n = case holds trivially. If we freeze «eZ\{0}, we find that 

If _2Trmx X - r- . TTjy 

— / q(x,y)e l dx — a n o + ) j a n k V 2 cos — — = ^ & n . ; v2sHi (3.9) 

fc=i j'=i 

Thus, the coefficients a n /c and 6 ra j are related to each other by a unitary transformation 



DC: 



a nk = E ^ b nj, (n e Z, fc > 0). (3.10) 

The entries of E can be computed explicitly: for j > 1 we have 

2v^/0V), fc = 0, j odd 



, Jo V2sm(Trjr))dri, k = 0\ j 4/ 

^fci = 4 ,i > = < , 9 7 , , fc>0. j-fcodd 

J 2sin(Trjr])cos(TTkr])dr], k>l) | (j fc J 71 " 

0, otherwise 

(3.11) 
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Keeping n € Z \ {0} frozen and dividing (|3 . 8|) by LH, we must show that 

Eoo (2-rm/L) 2 \a nk \ 2 < ^ ^oo (2-n-n/L) 2 \b nj | 2 . ^ -.x-^ 00 (-rrk/H) 2 \a nk \ 2 

fc=0 (27rn/L) a +(7rfc/_H") a - Ul = l (2im/ Ly 2 +(nj / H)* " t "l L/ 2 JJ Z^fc=l (27rn/L) ;i + (7rfc///) ;i ' 

This is accomplished via the following lemma using = 2\n\H/L and = 2H/L. □ 

Lemma 3.2. Suppose b e ^ 2 (N) and Zet a = Eb <E ^ 2 (No), w/iere _E maps the 
Fourier sine coefficients of a function to its Fourier cosine coefficients; see iS.ll]) 
above. Then for v > there holds 

oo 2 00 2 00 / 2 

fc=0 j = l J k=l 

with C\ = max (9, j^~ 2 ) and C2 — 9. If v > > 0, Ci = max (9, f^ 2 ) aZso works 

Proof. It suffices to show that p,12j) holds whenever b is a unit vector in i? 2 (N). 
The general case follows by re-scaling this result. We will split each sum into terms of 
low and high index and use different arguments to handle the two cases. Let fco > 0, 
jo > 1) k\ = fco + 1 and j\ = jo + 1. If we discard terms on the right hand side with 
j > ji and k < ko, we obtain a sufficient condition for (|3.12|) to hold. Also, on the 

left hand side, J^kLki T^Tk? \ a k? < Efclfei ^iwW? , so it suffices to show that 

fco 2 Jo 2 / 2 \ 00 j,2 

S CT |flt|2 * Cl S ^^ 2 + C2 - 1 - F ) E ^Trb^l 2 - ( 3 - 13 ) 

fc=0 7=1 J V 1 



fc=o i=i J v 17 fc=fc 

Next, we see that (13.13[) will hold if we can show that 

- ^TT!^ + ^Tk\ (1 a } ' (3 ' 14) 

where a 2 = ^fclo M"- ^ = Ef =1 l^| 2 , and 1 - a 2 = J2Z kl \ a k?- Note that (3 
here is not the inf-sup constant (3, but rather a measure of the relative weight of low 
frequency modes in comparison to high frequency modes in a sine series expansion. 
Solving (|3.14p for a 2 , we require 

gl l + , 2 /fc 2 l + , 2 /fc 2 

a ^c 2 i + Jo 2 /^ +1 c 2 ■ (3 - 15) 

Our goal is to show that for each v > there is a choice of jo > L &l > 1, C\ < 
max(9, (9/4)^~ 2 ) and C 2 < 9 such that ([3~15| and consequently (|3~T2]) holds for all 
unit vectors 6 e ^ 2 (N); (b determines a, a and (3). C\ and C 2 can then be increased 
if necessary to the values stated in the lemma without violating (|3.12j) . 

We now use the fact that a and b are unit vectors related by a known unitary 
transformation to obtain a bound on a in terms of (3. Let S, T, x, y, z be the 
sub-matrices and sub-vectors 



S = E(0: ko, l:j ), T = E(0:k , jf.oo), E(0:ko, :) = [S, T]. 
z = a(0:fc o ), x = b(l:j ), y = b(J 1 :oo), z = Sx + Ty. 



(3.16) 



We have a = \\z\\ and (3 = \\x\\ = \/l — ||y|| 2 - Since \\S\\ < 1 and \\y\\ < 1, the 
estimate ||z|| < \\Sx\\ + \\Ty\\ gives 

a<(3 + t, t=||T||<l. (3.17) 
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If t < 1, this can be used to derive a bound on a 2 of the form (|3 . 1 5[) . However, 
we can obtain a sharper estimate as follows. First, we compute the singular value 
decomposition S = UT,V* and rotate the rows of T by a unitary operator Q such 
that 



U*[S, T] 



V o" 




Q 





V o 















to 








tk Q 



(3.18) 



where a\ + 1\ 



1 for < k < ko- We assume here that jo > ko + 1; otherwise we will 
not be able to derive a sufficient condition for (|3.15p to hold, for if S has more rows 
than columns, we can produce a unit vector a = [z; 0] with S* z = so that b = E*a 
yields a — 1 and = 0. Next we define z — U*z, x = V*x, y = Q*y so that 



a 



En 2 ' 



Jo 



fc=0 j=l 

and, by Lemma lA. II below. 

1 . _ 



/5 2 = Ei%i 2 ' 



z k — a k x k+ i +t k y k+ i 

(3.19) 



\~z k \ 2 < 



\<TkXk+l\ 



1-tfc 

Hence, majorizing t k by ||T|| = t 



^-\t k y k +i\ 2 = (1 + t k )\x k+1 \ 2 + t k \y k+1 \ 2 . (3.20) 



J 2 . 
mm 



a 2 < (1 + t)(3 2 + t(l - f3 2 ) =(3 2 + t, 
Thus, (I3.15|) holds if we define C\ and Ci via ^1 — 



and summing over k, we obtain 

<=||T||<1. (3.21) 



c 2 



t and 



{C 2 1+JS/u' 



= 1: 



l - 1 



C 2 {v) = 



l + u 2 /k 2 
1 - t 



(3.22) 



Next we look for choices of k\ and jo that lead to a window of values of v over which 
C\ and C2 remain small. We need enough such windows to cover the positive real 
line v > 0. The trade-off is that choosing jo k\ makes t small but also makes one 
of the numerators in (|3.22p large. We consider 3 cases: 
• Case 1: (0 < v < 1/2). We set ki 



Jo 



1 so that t = y/l 



E oi 



.4352 and 



C 2 {v) = (1 + v 2 )/{\ -t)< (5/4)/(l - t) = 2.2133 < 9/4, 
d(v) - (1 + v- 2 )/(l - t) = C 2 v- 2 < (9/4K 2 , 



(0 < v < 1/2). (3.23) 



• Case 2: (1/2 < v < 100). We wrote a program to compute the singular value 
decomposition of S for all pairs of small integers k\ and jo satisfying k\ < jo < 4fci < 
400 to determine t = y/l — cr^ ln for each pair. We then choose a threshold Cthresh 
and find the values v min and f max such that Ci(^ m i n ) = Cthresh and C^^max) = 
Cthresh- We then discard all cases with f max < j/ min and sort the remaining intervals 
[fmin, i-'max] by their first entry. Finally, we discard all intervals for which v m in of the 
next interval is smaller than ^ ma x of the previous interval (to avoid redundancy) . The 
results with Cthresh = 4.9 and Cthresh = 8.9 are shown in Figure [3~T1 The method 
breaks down (i.e. there are gaps between some of the windows) for Cthresh < 5.83. 
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Table 3.1: Parameters used to 
construct C\{v) and C2{v) with 
Cthresh = 8.9. The correspond- 
ing table with Cthresh = 5.9 has 
32 lines corresponding to the 
smaller windows shown in Fig- 
ure 0Q3 



C^v) and C 2 (v) 



fcl 


Jo 


t 




^max 


1 


1 


.43524 


0.498 


2.007 


3 


3 


.57904 


1.810 


4.972 


6 


8 


.54892 


4.608 


10.42 


13 


17 


.58222 


10.31 


21.43 


25 


37 


.54766 


21.27 


43.49 


50 


76 


.54321 


43.41 


87.54 


99 


155 


.53535 


87.54 


175.3 



72 r 



max(9,(9/4)0 C lhresh = 8.9 



thresh 




Fig. 3.1: Plot of C\(y) and Ci(y) over the range 
0.2 < v < 600. Each criss-cross corresponds to a 
different window v m i n <v< f max in Table [3~T1 



• Case 3: (u > 100). We set k\ = \y/\/Z\, io = 3fci and bound t by the Frobenius 
norm: 



k() oo 

* 2 <imii = EEi^i 2 = ^E 

k=0 j=h 



8j,o. 



j, odd 



1G 



7T^ e — ' 7 

3=31 



■ , 1 , 8 , 
< — I —dx + — fc 

30-1 

4 



SEE 



j Sj-k,odd 



J (> 



< 



TT^Jo-l) 7T 2 

4 4 

1707T 2 + ^ 



-1 (X 2 -^ 
1, K+l 



1 2 ° K- 1 



7T^ ^ ' (? 2 - fc 2 ) 2 



(&c> (io - 1 = ji - 2) 



jo - i . io „ 

ft = — : > — = 3 



(3.24) 



3 1, „ 
- + - log 2 

8 2 S 



(.54298) 2 





100 


(fcl > 


L VsJ 



fc fcl 

57, io > 171 ) . 



Here we represent sums of decreasing functions sampled at even or odd integers by 
staircases of width two and half the height of the function at the right endpoint. Each 
choice of k\ and io will cover the range V3fci < v < \/3(fci + 1); over this range, we 
have 



v / fci + 1 
fci " 



fci 



—\ < SV3, t < (k < (3)^ = V3 (3.25) 



fci + 1 



57 



k\ J \v 



V3 



and we learn that C\ and C 2 are bounded by ^^s^gg = 8-985. 

Thus, for all j/ > we have Ci < max(9, (9/4)j/~ 2 ) and C2 < 9, as claimed. □ 

Corollary 3.3. For all q e L 2 # (R), ||g|| 2 _ 1 < ] | [ 2 _ 1 + fr||0 v ?||!. r 
Proof. Arguing as in (|3.4[) - (|3.6[) . it is readily shown that 



-x = E 



ZxN 



E 



LH\b 0j \ 2 



The first term on the right hand side is simply while the second satisfies 



E 



LH\b 0j \ 2 < H 



^ (nj/H) 2 ~ n 2 ^ 



J2 LH \ b < 



H 2 - 



Y^LH^ 2 <^\\d y q\\ 2 _ v (3.26) 



i=i 
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where the middle equality follows from the fact that aoo = 0. □ 

-2 



Corollary 3.4. Suppose q G L\(R) and £ G L°°(T). Then for any aspect ratio 



H/L. we have 



\\d y (Cq)\\ 2 ^ < C 2 M 2 (WdrfWU + Rgll 2 -!) 



(3.27) 



where C2 = 9 and Af = ||£||oo- 

Proof. Since C does not depend on y, the Fourier coefficients of q — (^q are related 
to the those of q via column- by- column convolution with the Fourier coefficients of £: 

a-nk = ^ (n-"i a mk, Kj = ^ Cn-m&mj) (n € Z, fe > 0, 7 > 0). (3.28) 



m£2 



m£2 



Since multiplication by £ is bounded in L 2 (T) by M, convolution with £ is bounded 
in ^ 2 (Z) by Af. Thus, by ([3^4]) . we have 



\d v (Cq)\\ 



< 



^EE 



5 Ilfc | 2 < M 2 



(3.29) 



fe>0 nGi 



fc>0 n£2 



The key point is that entries a n k with k = are absent from the right hand side. The 
quantity in parentheses may be written as A\ + A2 just as in (|3.4[) . but omitting the 
k = terms from A%. Thus, it suffices to show that (|3.12D holds with C\ replaced by 
Ci if the k = term is omitted from the sum on the left. For v > 1, the result has 
already been proved without omitting this term. But for v < 1, we see that C\ = 
and C2 = 2 suffice (since ^ 2 < k 2 for fc > 1). Thus C\ — C2 — 9 works for all ^ > 0, 
as claimed. □ 

4. Curved boundaries. We now perform a change of variables to transfer the 
result of Theorem 13.11 from the rectangle R to a domain SI bounded on one side by a 
periodic, Lipschitz continuous function h G C°' 1 (T): 



7/ 





h 






F 






0. 


Q, x, y 



H 



y, d(, dr) 



H 
H 

h(x] 



dx dy, 



d d 77 d d H d d d y d d h d 

dx d£ h x dr] 1 dy h dr)' <9£ dx h x By' di] H By ' 

(4.1) 

The main challenges involve avoiding lower order terms that have to be dealt with 
using Rellich's compactness theorem, balancing the sources of error to avoid excessive 
overestimation of the constants in the error bounds, and dealing with various subtleties 
of the dual space if -1 (SI) such as the fact that if p £ L 2 (S!) and C £ L°°(S7) then 
||Cp||-i need not be smaller than ||C||oo||p||-i- For clarity, we postpone the case that 
h is only Lipschitz continuous to Section [5] and begin with the simplifying assumption 
h G C X,1 (T). The aspect ratio of the rectangle R plays an essential role in the Lipschitz 
case but only a minor role (improving our estimate of j3) here. 

Theorem 4.1. Suppose h <E C 1,:L (T) and < h < h(x) < hi for < x < L. 
Then for every p £ i^(Sl) we have 



lo.fi 



-i,fi' 



0- 



-(1 

4 V 



m 2 ) (h 



h 



1/2 



h h J 
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where M* = ^(\\h4lS^\L)- 

Remark 4.2. The quantity \hh xx arises naturally in the study of Reynolds' 
lubrication approximation and its higher order corrections on a periodic domain 12]. 

Remark 4.3. In many practical applications, the aspect ratio L/ho is large while 
M<1 and hi /ho ~ 1; in this regime, (14. 2| shows that j3 scales linearly with L/ho. 
If the geometry has a narrow gap so that hi /ho ^> 1, we learn that depends 
on the gap size as h 3 ^ 2 . This dependence is shown to be optimal in Example 14.51 
below. We do not know if the quadratic dependence on M is optimal; it seems to be 
an unavoidable artifact of changing variables to a rectangular geometry. 

Proof of Theorem \4-l\ The coordinate transformation (x,y) = F(J[,rf) defined in 
(|4.ip provides a one-to-one correspondence between functions p G L 2 (il), u S Hq(H) 
and their counterparts p — p o F € L 2 (R), u = u o F 6 Hq{R): 

V) = P (t, ^) , fi(£, *?) = « ^) • (4.3) 

F does not map i|t(fi) to L^(i?); however, the norm of p G L^(fi) does not decrease 
if we add a constant to p to enforce § n h~ 1 p dA = instead of J n pdA = 0. By 
Theorem 13.11 this new p satisfies 



I 2 < 
lo.fi — 



/^i\ 1/2 " ftl IMI 2 ^ n hl \\x ~\\ 2 ,^< ftl ll«~ll 2 ,n 
(XJ P o,fi = ^IWIo,K< ^^11^11-^ + Ca^HM-x,*. (4-4) 



where Ci = max(16, -^p-) and C 2 = 9. But since the right hand side does not change 
when a constant is added to p, the original p also satisfies this equation (dropping the 
intermediate inequalities). We can relate the action of V^p on u to that of V x p on u: 



R = (jjp, - \ h * d y) M ^ = H (9xP, h 1 u) n + H (d y p, h x j^uj^ . 



(4.5) 

(d v p,v) R = (^P'-Ji v v/ ={ 9 vP' v )ni (4- 6 ) 

where we used d x (h~ r ) + d y (yh~ 2 h x ) — in (|4.5|) . If we had not introduced the 
factor of h~ x l 2 in (14. 4|) . this cancellation would not have occurred and the proof 
would become much more complicated; see Remark 14.41 below. It will be shown in 
Lemmas IA.31 IA.4I and IA.5I that 

H\\h-^u\\ a ^<C*\\u\\ a ^ C 3 2 =max^£, (l + 3M 2 )|J), (4.7) 
< C* 4 ||u|| ait , C\ = max [%M 2 — , (2M 2 + 6 A/ 4 ) J, (4.8) 



H 



h x - 2 u 



HLfi< C 5||*IU> Og=max^, (l + 2Af 2 )£ j. (4.9) 

If /i only belongs to C 0,1 (T), then (|4.8|) docs not hold and we have to replace the last 
term in (|4. 5[) by (d y (h x p), Hyh~ 2 u) w which requires a more difficult analysis; see 
Section [5] below. Combining (|4.5j) - (|4.9p . we obtain 



\{d ( p,u) R \ < (C 3 \\d x p\\_ in + C i \\d y p\\_ in j\ 
\(dr,P,v) R \ < c 5 ||a y p||_ 1 n ||«|| fl . 
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It follows that 

\\d ( p\\ 2 _ 1>R < 3C|||5xP||!. lin + \cl\\dyP\\\sv 
which, together with l|4.4j) . gives 



/'"'(lIMI-x.n+IIMI-x.n). ^ 
Next, we choose H = ho so that 

3CiC 3 2 < 9(1 + Af 2 )Ci, |CaC2 < (l2Af 2 



hi 
H 



max SdCl, -dC 



9Af 4 )C!, C 2 Cf < 18 



c 2 c 2 



18Af 2 



Finally, we observe that ^ < ~ max ^16, j^J regardless of whether ^ > 4. As a 
result, C 2 Cf < (8 + 2A/ 2 )^ max (l6, j|) 



and 



h 9 / L 2 h? 

/T 2 < — - max {9(1 + Af 2 ), 8 + 14Af 2 + 9M 4 } — max 16 —,-1 
ho 16 \ h h 



(4.12) 



which yields f|4. 2|) when we majorize the terms in braces by 9(1 + M 2 ) 2 . □ 

Remark 4.4. One might hope to improve (|4.2p by working directly with \\p\\o 
in l|4.4[) instead of via ||/i _1 / 2 p|| . The main difference is that (|4.5|) acquires a lower 
order term 

(0*(fc 1/a p), u) R = H(d xP , h~ 1 ' 2 u) n + H(d yP , yh x h' 3 / 2 u) n + jfa h x h- 3 / 2 u) n 

that would normally be dealt with by invoking a compactness argument to bound 
||p||_i.n by a constant times ||Vp||_ i,n- This is not acceptable in the current calcula- 
tion as this constant depends on f2, and hence h. It is possible to bound x,n in 
terms of ||p||_i,ij and then use Corollarv l3.3l But the final step of bounding || VgpH-i.ij 
by llV^pll-i.n brings us back to the proof given above. The following example shows 
that the power of h 3 ^ 2 in the formula (|4.2[) for fl^ 1 is the best possible. 

Example 4.5. Suppose < ho < 1 and consider a periodic function h(x) that 
transitions smoothly and symmetrically between ho for x € [3/8, 1/2] U [7/8, 1] and 1 
for x G [1/8, 1/4] U [5/8,3/4]. Let fii, f2 2 , and be the regions under the curve 
h with x G [0,3/8], x G [3/8,1/2], x G [1/2,7/8] and x G [7/8,1], respectively. Let 
p(x,y) be the continuous, piecewise linear function that equals —1 on fii, 1 on ^3, 
and satisfies p x = ±16, p y = on 17 2 and O4. 



1 

h) 








\p x = \Qf \p, 


/ 




\_j V 




p 




1 Q 2 p = 1 0. 







Then for any u £ H q (£1) , we have I (d y p, u) I = and 



\(d x p,u)\ < 



16\u(x, y)\ dA < 16y / area(fi 2 U Q4) ||u||o,n 2 un 4 



< 



o 2 ur2 4 
1/2 



u y \\oM 2 uni < Vshl \\u\\ a ,n, 



(4.13) 
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where we used the Cauchy-Schwarz and Poincare-Friedrichs inequalities; see Lemma 
|A~21 Thus ||Vp||_ in < V&hl /2 while ||p|| ,n > 1/2, showing that ft- 1 in (O must 

be at least (2\/8) h Q , i.e. the power h 3 ^ 2 is optimal. 

5. Lipschitz boundaries. In this section we show how to modify the proof of 
Theorem 14. II to handle the case that h only belongs to C 0,1 (T). The main difference 
is that yh~ 2 h x ii no longer belongs to Hq (fi) in (|4.5p . so a different strategy is required 
to deal with the term (d y p, Hyh~ 2 h x u)^. The idea is to show that when h~ 2 h x is 
grouped with p, this term can be made small in comparison to the other two terms 
in (|4.5[) by choosing the aspect ratio of the rectangle R small enough. The loss of a 
power of h^ 2 in the estimate of when M is not small is discussed in Remark 1 5. 3 1 
below. 

Theorem 5.1. Suppose h £ C 0,1 ^) and < h < h(x) < h\ for < x < L. 
Then for every p £ Lg(Sl) we have 

IMU * n\Vp\Um = 2 max (4, ^,8^m) max(l,8M)^ 

(5.1) 

w/iere M = ||/i x ||oo- 

Proof As before, (|4.4[) holds for all p e L#(fi): 

\\p\\l,n^ c ^\\ d M\ R + c ^\M\-i,w °i = Ye max ( 16; §2 ) . C2 = 9 - 

(5-2) 

We now transform the problematic term in (|4. 5|) back to the £, ?y coordinate system: 

(f,u) R ~ (9i,u) R := (d^p,u) R - (d v (h~ 1 h x p),i]u) R = (d x p,Hh~ 1 u) n , (5.3) 
(ff: 5 >fl : = {d v P,v) R = (d y p,v) n . (5.4) 

So we can bound ||f>||o,n hi terms of / and g and we can bound (/ — g%) and g in terms 
of || V x p|j_i^n; thus, we need a bridge from / to (/ — g{) and g. By Corollary [33] and 
Lemma IA.6| 

\{gi,u) R \< \\d n {h- x h x p)\\- l<R \\m\\ a ,R< (sM/^iiVepll-i,*) Q#INkfl) , 

||5iH 2 -i, fl < ^(ll/lP-i,* + Ibll-iJ , = 4^M. (5.5) 

As a result, ||/|| 2 < 2||/ - , 9 i|| 2 + 26 2 {\\f\\ 2 + || 5 || 2 ), which implies 

||/|| 2 <4||/- 3l || 2 +40 2 || 5 || 2 , (tf 2 <l/4). (5.6) 
Equation (|5.2[) now becomes 

Ibllo.n < 4Ci^||/ - ffl || 2 _ l!R + (40 2 d + C7 2 )^|| ff || 2 _ l!R , (M < fto/8). (5.7) 
From (|5.3p and ([5T4]) we see that 

|(/-ffi,«)|< ||^||_ 1) n||iI/i- 1 w|U,n 1 |{fl,5)|< Rp||-i,n||«||«,n. (5.8) 
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By Lemmas IA.3I and IA.5I below, we then have 

11/ - 9x\\-l,R < C 3 \\d xP \\- h n, Cl = max ( ~, (1 + 16Af 2 )^ ) , (.,.«.) 



8 ho ,y ' h\ 

MI-i./.-l C:\\i) ul >\U.>.i. Cl = max (£^,(1 + 9M 2 )^V (5.10) 



It follows from (15.71) that 



' 2 <r 2 (||MI-i^ + IIMI-iJ' r 2 = ^-max(4C 1 C 3 2 , (40 2 d + C 2 )C 2 ; 



|Pho.n 



Finally, we choose H = min (h$, g^y/io) so that if M > 1/8 we have 

hi „ / 9 4 \ fel fel ftl „, / ft? 9 . 0>ftl\ 9 ft? 

4— Cl <max -, 7 + 1) — <5 — , — C? < max 72 -| Af 2 , 1 + 9M 2 — <73M 2 ^- 

H 6 ~ V 2 64Af 2 J ho ~ ho H \ h o h J ~ ft 2 , 

and if Af < 1/8 we have 

hi 9 /9 ,\ hi fti fti , /9/i? , 9^1 \ 73 ft? 

4— Cl < max -,4 + 64A/ 2 — < 5—, — C? < max --i, (1 + 9M 2 )— < L 

H \2 J ho ~ ho H 5 ~ V8ftg / ft y~64ft 2 

Moreover, Ci = max (9, 36^M 2 ) and 46> 2 Ci + C* 2 < 2 max (9, 36^M 2 ) re- 

gardless of whether M < 1/8. Combining these results, we obtain 

4C 1 |c 2 <5max(9,l|,36|M 2 )^, 



)^C 2 < I max (9,36-gA/ 2 ) max(l, 64M 2 )|. 



(5.11) 



Formula (|5.1[) for /3 _1 follows by taking the square root of the maximum of these 
expressions after increasing the constants and consolidating terms. □ 

Remark 5.2. Inequality (|5.5j) is the key to this proof. For fixed it, both (/, u) 
and (gi,u) in (|5.3|) scale like 77 while in (|5.4p is independent of H. Because of 

the way ||u|| a ,fl depends on H, it follows that if R\ = T x Hi, R2 = T x H%, and 
Hi < H 2 , then 

\\9l\\\ Rl < ^IM| 2 -l,H 2 , ll/ll 2 -!,^ < ^ll/l| 2 - 1)A , llffll 2 -^ < ^\\9\\\ R2 - 



Thus, ||<7i|| 2 and # 2 ||g|| 2 are both 0(H) quantities and the surprising aspect of (|5.5|l 
is that the 0(H 3 ) term 6> 2 ||/|| 2 is sufficient to help 9 2 \\g\\ 2 bound ||5i|| 2 . 

Remark 5.3. We believe the optimal bound in the Lipschitz case should scale like 
/3 _1 ~ h Q 3 ^ 2 , just as in the C 1 * 1 case; however, proving this would require eliminating 
(or at least finding a better bound for) the cross term 4# 2 Ci^||g|| 2 _ 1 R in (|5.7|) . As it 
stands, 6> 2 Ci and i? _1 ||(7|| 2 each contribute a factor of h^ 2 to this cross term due to 
the requirement HM < ho/8, which yields j3~ 2 ~ h^ . We suspect that the functions 
p that require Ci to diverge as H — > are distinct from the functions p for which 
11/ — .9i II ^ 11/11 m l|5.6p . but we have not found a way to make this idea rigorous. 

Appendix A. Useful Lemmas. In this section we gather several results that 
are either elementary but used frequently in our proofs or are tedious and distract 
from the main argument. 

Lemma A.l. Suppose 71, . . . , y n are positive real numbers such that jj 1 < 1. 
Then \wi + --- + w n \ 2 < £™ "fj\wj\ 2 , for all w G C". 
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Proof. This is a consequence of the Cauchy-Schwarz inequality: 



-l/2\ / 1/2 
1 '.Tj 



< 



(A.l) 



Lemma A. 2. (Poincare-Friedrichs inequality) . If f2 has the geometry of Fig- 
ure \2.1\ with h S C 0,1 (T) and i/ R is the x-periodic rectangle of width L and height 
H , then 



h 2 

Mlo,n < -^-\\ u y\\ln, 



H 



In < -Fll%11o.il. (u G H^Q), « G ^(fl) 



(A.2) 

T/ie former inequality also works over the subregion VL2 U O4 m Example \4-5\ with h\ 
replaced by the maximum height of that subregion, namely ho . 

Proof. The latter inequality follows by expanding u = E d n j y/2 exp( 27r ^ na; ) sin 
and comparing the formulas for |M|o fl and ||u y ||o fl . If u G (7^(0), the former in- 
equality follows by integrating 



\u( Xl y)\ 2 < \J v u y (x,y')dy'\ 2 < y J Q h/2 \u y (x, y')\ 2 dy' , (0 < y < \h{x)) 

\u(x, y)\ 2 < | u y (x, y>) dyf <{h-y) jj /a \u y (x, y')\ 2 dy', {\h{x) < y < h(x)) 

over the lower and upper halves of Q, respectively, and combining the results. The 
result for u E Hq(&) then follows by a standard density argument. □ 

Lemma A. 3. Suppose h e C° 4 (T), R = T x H and u G H^n). Then 



\Hh-i U \\l n <Cl\\u\\l R , 



(A.3) 



where u expresses u in the £, rj coordinate system of R and Cf is given by 7[ ) 

/op . 

Proof. Using the change of variables formulas (|4.ip and (|4.3p , we obtain 



\d x (Hh~ l u)\( > Qn = J H 2 (-h- 2 h x u + h- x ut: - h~ ir ^h x u n 



2 h 
H 



d£ drj 



< 7l ifV 3 ^ 2 INl6U + 72Hho 1 ||uc||S,R + 73^0 3 Af a ||fi,||g lR , 



3 1,-3 ji.j-211 ~ ||2 



(A.4) 
(A.5) 



where M = \\h x \\cx>i ^0 = niino< x <L ft(ic), and 7l 1 + 7 2 1 + 7 3 1 < 1. Combining these 
and using the Poincare-Friedrichs inequality (with 9 instead of 7r 2 ), we find that 



\Hh 



-i„.ll 2 



H 



7i 



„,„ < max ( 72 -, (1 + ( 73 + J± j AT' J ^ ) ||n||^, (A.6) 



which yields (gTTJ) with 7 = (9, 3, 2) and (EU) with 7 = (36, 9/8, 12). □ 

Lemma A. 4. Suppose h e C M (T), R = T x H and u G ™en 

||iTJrV»**|lln < C 4 2 ||u||^, C 2 = max fsM 2 ^, (2Af 2 + 6Af 4 )|J ) , (A.7) 



h 3 
"0 
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whe n J^ = mBx(j\h m \\l^hh am \\ go ). 

Proof. Using the change of variables formulas (|4.ip and (|4.3[) as well as the 
Poincare-Friedrichs inequality, we obtain 



0x( 



Hyh x u \ 2 
h 2 J o,Q 



u + 2 _^ftft xa j u +_^ , 



2 h 

— d^dr] 



< 4(7i +72) 
HM 2 



HM 



HM 2 



h 3 n«llo 1 H + T3-^-||^llo 1 ii + 74 

H 3 M 4 



"■0 



73 _ /^ _ " u? " '- R + V 74 + ^2 (71 + 72) j fe3 IKHo.fl ( A 



r Hyh x u \ z = I ( 
\ h 2 J a si In V 



i/Zij. _ Hr)h x \ 2 h 



h 2 



-u + 



h 2 



H 



d£ drj 



(A.9) 



"0 



ft 



H 3 M 2 , 



.3 ll^nllo.-R' 
"■0 



'9 ip 
i2' 2/ 



where Jjj 7," 1 < 1 and #1 < !• Now we set 7 = (§tt 2 , |tt 2 , 8, 3) and 6 

to obtain (|AT7|) . □ 

Lemma A. 5. Suppose h g C - 1 ^), R = T x H and v g Then ||u|| a ,n < 

C5||w|| ,fl wfi C5 ffwen &y gT3|) or jj.ifl)) . 

Proof. Let M = ||/i x ||oo. F° r & n y Tij 72 satisfying 7f 1 + 72 -1 < 1, we have 



\\d x v 



2 

lo.n 



_ \ 2 



„) ^dfd»?<71^IKIIo,B + 72-^-|Kllo,fl. (A.IO) 
\\dyv\\l u = j v v f^dr,< £lNlo,n- (A.ll) 



It follows that ||«|| a ,n < C 5 \\v\\ a ,R with C 5 = max (71^, (1 + j 2 M 2 )^). We obtain 
(|4~9]) using 7 = (2, 2) and dSHUJ) using 7 = (9/8,9). □ 
Lemma A. 6. On the ^-periodic rectangle R, 



16 



hcC<-# 2 ||tC, (ie^)), 



(A.12) 



Proof. Using the Poincare-Friedrichs inequality, we find that 
\\ddm)\\l, R <H 2 \\u 6 \\l R , \\d v ( V u)\\l R < (|-+72) J ff 2 |KlliU, (A.13) 

provided 7-f 1 -I- 7J < 1. Choosing 7 = (4,4/3), the result follows. □ 

Appendix B. The Poincare-Friedrichs inequality on ffi(O). In this section 
we present a simple proof of the Poincare-Friedrichs inequality for H 1 functions with 
zero mean. Our proof does not rely on Rcllich's compactness theorem, but does 
require the boundary of O to be the graph of a Lipschitz continuous function h; see 
Figure [2~T1 above. The main difference between the estimates 

IWko < A-||Vp|| ,n, (P eH U Q ))> lbllo,n<j8- 1 ||Vp||_i, nj (p G L%(Q)) 

(B.l) 

proved below and in Theorems l4.1l and l5.1l above is that K ~ h Q 172 while fi^ 1 ~ h Q 3 ^ 2 ; 
(we were only able to prove (3~ x ~ h^ 2 in the Lipschitz case). A narrow gap causes K 
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to grow because a large gradient of p in the gap region can lead to a large change in p 
across the gap with relatively little cost (in terms of ||Vp||o,o) due to the small area 
of the gap region. The effect on is more severe than on K because, in addition 
to the small area of the gap region, the test functions (u, v) that Vp acts on belong 
to Hq(H) 2 , i.e. they are zero on Tq and Ti. These boundary conditions cause u and v 
to be small in the gap region, which reduces their ability to penalize large gradients 
of p there. This was illustrated in Example [43] above. 

To keep the equations dimensionally correct, we define the norm on H^(fl) to be 

IWIi.n = L- 2 \\p\\ln + \\p\\l,n = J ^ + W 2 + \Py\ 2 dxdy, (B.2) 

i.e. we use L as a length scale to compare ||p||o to ||p|j Q = ||Vp||o- 

Theorem B.l. Suppose h E C°^(T) and < ho < h(x) < hi for < x < L. 
Then for every p E H^(fl), we have 



i-'IMkn < C||Vp||„,n, 



C = 



l + M 
2tt 



max 1,2 



L 



lh[ 
h ' 



(B.3) 



where M = \\h x \\oo- The constant K in \B. 1\) is given by K = (1 + C 2 ) 1 / 2 . 
Proof. On the ^-periodic rectangle R — Tx (0,H), the expansion 



irkn\ 2 ""£ 



H 



P(€> V) = E ( a "° + E a nkV2cos 

riGZ ^ fc=l 

can be differentiated term by term and we have 



EH\R)) (B.4) 



/2irn\ 2 / nK 



rk\ 2 



v L , « „ ; ■ (B.5) 



H 



wpwIr = E LH \ a n k \ 2 , \\vp\\Ir = J2 lh 

n,k n,k 

Assuming p £ i.e. aoo = 0, we learn that 

\\p\\l R < L 2 C 2 \\Vp\\t, R , L 2 & = max {(A) 2 , (^} 



\a n k\ 



(B.6) 



Now we transfer this result to VL by the change of variables (|4.1|) and (|4.3[) . To avoid 
Rellich's theorem, we estimate 



This inequality holds for all p such that p E H^(R). Arguing as in (|4.4| . we find that 
if we drop the intermediate inequalities, (|B.7[) also holds for p E H^(fl). Next, we 
bound ||Vj5||o._r, in terms of ||Vp||o.j?: 



\\Pt\\o,R = 



\\Pv\ 



Px + ^hxPyj — dxdy < 71 — \\p x \\o,n + I2M ^||p y || Q) n 

2 f ( h \ 2 H , , hi, 
> r = J u {h P «) -h dxdy ^H^ 



2 

o.n- 



where r y 1 1 + 7 2 1 = 1 and M 

L~ 2 \\ P \\ln<C 2 \\Vp\\ln, <' 



211. ,12 . ,-11^,112 ■ -C^mzxU 1 -!- '" ' - ^ 

V h Q 



H hi ir0 H 

ff max l 7l ^'77 + 72M h a 



(B.8) 
(B.9) 

(B.10) 
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Next, we choose H = \/h hi and minimize max (71, 1 + 72 M 2 ) over all choices of 7^ 
such that 7-f 1 + 73 = 1- The result is 

71 = l + 72 Af 2 - ^(y/AP + 4 + M) 2 < (1 + M) 2 , (B.ll) 

which yields C 2 = ^ max (l,4^i) |i(l + M) 2 as claimed. □ 

Remark B.2. Example 14.51 shows that the scaling C ~ /i is optimal: for 
that function p, we have 




(B.12) 



which shows that C in (|B.3|) is at least j§\J j^- We do not know if the linear depen- 
dence of C on M is optimal — it seems to be an unavoidable artifact of changing 
variables to a rectangular geometry. 
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